Non-coalescence of sessile drops from different but miscible liquids: 
Hydrodynamic analysis of the twin drop contour as self stabilizing, traveling wave 
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Capillarity always favors drop fusion. Nevertheless sessile drops from different but completely mis- 
cible liquids often do not fuse instantaneously upon contact. Rather, intermediate non- coalescence 
is observed. Two separate drop bodies, connected by a thin liquid neck move over the substrate. 
Supported by new experimental data a thin film hydrodynamic analysis of this state is presented. 
Presumably advective and diffusive volume fluxes in the neck region establish a localized and tem- 
porarily stable surface tension gradient. This induces a local surface (Marangoni) flow that stabilizes 
a traveling wave i.e., the observed moving twin drop configuration. The theoretical predictions are 
in excellent agreement with the experimental findings. 

PACS numbers: 68.03.Cd, 68.03.Kn, 68.08.Bc, 83.50.Xa, 83.60.Uv 



Introduction. — Basic physics tells us that two sessile 
drops will coalesce as soon as they touch each other be- 
cause of the reduced interfacial area/energy for a single 
drop. For drops with identical liquids the instantaneous 
coalescence after contact has been studied [HQ. Yet, lit- 
tle is known about the fusion of drops with different, but 
completely miscible liquids. Recent studies show that it 
can be fundamentally different |3|— ISI] - The drop fusion 
can be delayed for a long time. After contact, the drop 
bodies remain separated in a temporary state of non- 
coalescence, connected only through a thin liquid bridge. 
This "twin drop configuration" moves over the surface 
(Fig. 1). Presumably a Marangoni flow resulting from 
the surface energy difference of the two liquids causes 
the non-coalescence [3|, |4[ . An adequate hydrodynamic 
description of the non-coalescence state has not yet been 
presented. Especially, it has not yet been explained how 
a temporarily stationary Marangoni flow is established, 
how this flow stabilizes the non-coalescence, and how this 
moves the twin drop system. 

Understanding the phenomenon is of general interest 
because it is counterintuitive, seemingly ignoring capil- 
larity. The interaction of drops of different composition is 
relevant e.g. to droplet-based microfluidics @, Q and es- 
pecially to Marangoni flows in thin films that are widely 
applied for surface cleaning ("Marangoni drying", 
Here we present new experimental data and the first an- 
alytical hydrodynamic description. We show that the 
moving twin drop configuration is a special case of tem- 
porarily stationary, self-stabilizing thin film waves that 
are driven by stationary local surface tension gradients. 

Experimental Results. — Figure [1] shows A) schemat- 
ically the experimental setup; B) the neck in de- 
tail: capillarity fills the neck from both sides whereas 
Marangoni sucks out liquid asymmetrically, thus mov- 
ing the drops with vn and preventing coalescence [lij : 
C) non-coalescence; D) instantaneous coalescence. Drop 
1 is deposited with a syringe from the top, drop 2 is 
formed by pumping liquid 2 through a hole in the sub- 
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FIG. 1: A) Schematics of the experimental setup/procedure; 
B) Flow scheme of the neck region; C) non-coalescent twin 
drop movement (different, miscible liquids); D) instantaneous 
coalescence (identical liquids). 



strate. More experimental details are found in [4]. With 
identical liquids, both drops fuse instantaneously upon 
contact (D). With sufficiently different surface tensions 
(C), the drops do not fuse [15j. The main drop bodies 
(h"(x) < 0) stay separated in a temporary state of non- 
coalescence, connected by a thin liquid neck (h"(x) ^> 0). 
Through the neck small amounts of liquid continuously 
flow from drop 2 (lower surface tension, 72) to drop 1 
(higher tension, 71). The flow slowly reduces the surface 
energy difference A7 = 71 — 72. Eventually the drops 
merge (delayed up to minutes |4|). The non-coalescing 
twin drops move over the substrate with constant veloc- 
ity, almost independent from A7 (Fig. [2]). 

The flow patterns within moving twin drops were in- 
vestigated by imaging dispersed fluorescent polystyrene 
microspheres (Duke Scientific, d = 1 |im, mass fraction 
w 2.4 • 10 -7 ). Figure [3] shows their traces in the neck 
region [14]. In the substrate frame all microspheres move 
from left to right. At /near the liquid/air interface of 
drop 2 the spheres move with « ^v N (positions 1 and 2). 
Close to the neck they touch the substrate (d « neck 
height h) and slow down (position 3). On the other side 
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FIG. 2: Neck displacement velocity Vat as function of the sur- 
face energy difference A7 / 72 (the insert demonstrates the con- 
stancy of v N ). Dashed lines: analytically calculated vn (see 
also Fig. [5}, in quantitative agreement with experimental re- 
sults. 



of the neck, in drop 1, those spheres that where at rest 
close to the substrate (position 4) get uplifted by the 
liquid flow through the neck. Spheres which reach the 
liquid/air interface of drop 1 can get accelerated up to 
~ 2v N (position 5). While moving away from the neck 
region, their speed slows down to « ^v N (position 6). As 
discussed below, the speed component in addition to 1^ 
is caused by the local Marangoni flow. 

Hydrodynamic model. — The steady state twin drop 
movement is analyzed by assuming a viscous, laminar 
Newtonian flow (no slip, no gravity, lubrication approxi- 
mation). The two drops are approximated by cross sec- 
tions through two infinitely long cylinders. They are 
connected by a thin neck and move slowly with con- 
stant speed and stationary shape over a planar substrate. 
They consist of miscible liquids with different surface ten- 
sions. Liquid with the lower surface tension continuously 
flows through the neck into drop 1. Diffusive- ad vective 
liquid mixing establishes a local stationary composition 
profile because the volume exchange per time is small 
compared to the overall drop volumes. The composition 
profile causes locally a surface tension gradient on drop 1 
close to the neck (for a detailed qualitative description 
how the Marangoni flow resulting from this gradient sta- 
bilizes stationary non-coalescence see [14]). 

The (stationary) surface topology h(x) of a liquid film 
with surface tension j(x), moving on a solid support with 
capillary number Ca = 77 v N /j is described by the balance 
between changes of the surface curvature, viscous dissi- 
pation and surface energy gradients, respectively [9]: 



3Ca/h 2 -377(2/17), 



(1) 



Without Marangoni term (7' = 0), Eq.Q] cannot describe 
a stationary twin-drop profile because h!" > allows only 
one inflection point (one drop connected to a neck) . With 
7' ^ however, a second inflection point respectively a 
second drop can exist because h!" may change its sign. 
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FIG. 3: Flow velocities (substrate frame) of individual mi- 
crospheres as function of their distance to the neck (moving 
twin drops, open circles: upstream, drop 2; closed circles: 
downstream, drop 1). Dashed line: measured neck velocity 
vn- Speeds exceeding ^vn indicate a (local) surface tension 
gradient. 

The velocity field u{z) within the film is 0|: 

u = z/n(j' - (z/2 - h)jti"). (2) 
Eqs. [T]and [2] yield the speed at the free surface (z = h): 



3^/2 + h 7 7 (4/y). 



(3) 



The Marangoni component increases u s : it sucks liquid 
out of the neck region and thus favors non-coalescence. 

This describes the experimental results (Fig. [3j). For 
drop 2, the maximum (surface-) velocity is « ~v N (= u s 
with 7' = 0). There is no surface tension gradient on 
drop 2. At the surface of drop 1, close to the neck, how- 
ever u s > ^v N : there, 7' 7^ 0, obviously caused by liquid 2 
flowing through the neck into drop 1. The data show that 
7' 7^ is localized next to the neck 0.2 mm, compared 
to drop sizes > mm). 

Localized surface tension gradient. — Marangoni forces 
cause the flow through the neck: they "pull" liquid 2 
as a thin film onto drop 1 where it mixes with liquid 1. 
This process is approximated by a layer of liquid 2 with 
the thickness of the neck, h N ,spreading on the surface of 
drop 1 with a velocity 1^ [16|. As liquid 2 flows onto 
drop 1 it is increasingly diluted by diffusion (diffusion 
constant D). Compared to the influx of liquid 2, drop 1 is 
a large reservoir of liquid 1. The dilution is a continuous 
and approximately stationary process. 

With negligible diffusion in flow-direction, this is de- 
scribed by the stationary advect ion- diffusion equation. 
In the (moving) contour frame the advection velocity is 
v N /2: 



d 2 c 



(4) 



dz 2 2 <9x' 

with local compositions c G [0,1] from pure liq- 
uid 1 to pure liquid 2, respectively. c \ x=0 = 1 f° r 
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FIG. 4: Moving drop contours with different shape param- 
eters S2- Dashed black lines: analytical solutions for homo- 
geneous surface tension (single drop). Solid red lines: non- 
coalescing twin drop contour (numerical solution) for a local- 
ized surface tension gradient 7' (top panel), and the corre- 
sponding velocity field (middle). 



< z < 2h N and elsewhere (2 ft, assures dc/dz\ z _ h = 
i.e., non- volatile liquids). j(x) shall depend locally on 
c(x,z) with j(x) = 72 + Ayc(x,ft). With this the solu- 
tion to Eq. |H yields: 
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Bo = v N h N /D is the Bodenstein number, a Peclet Num- 
ber with advection and diffusion orthogonal. 

i max is at x max Boh N /12 with Ajo fwhm ~ 

0.225Boh N (Fig. |4j). With typical experimental param- 
eters (v N ~ lmm/s, h N ~ 10 |im, D « 10~ 10 m 2 /s) 
7' is ^ only close to the neck (x max « 0.08 mm, 
Axfwhm ~ 0.2 mm). This and typical experimental val- 
ues of A7 ;$ 0.1 72 justifies to assume 7(2;) « 72 whenever 
j(x) is used explicitly. 

Numerically determined twin drop contour. — Eq. [T] to- 
gether with Eq. [5] is solved numerically. The result 
(Fig. 31 solid/red lines) is remarkably similar to the ex- 
perimental findings (explicitly shown in 0, 0). It in- 
deed reveals a moving twin drop configuration driven by 
the localized surface tension gradient described by Eq. [5] 
The contour for x < (drop 2), including in particular 
^"la^o * s described by Eq. [T] with 7' = 0. The main 
body of drop 1 also has a contour very similar to that of 
drop 2 i.e., with 7' = 0. This and the limitation of the 
gradient to the neck section is the key to the following 
analytical description which aims at predicting v N from 
only measurable quantities. 

Analytical approach. — With the localization of 7' 7^ 
at the neck region of drop 1 the remaining bodies of 
the two drops can be described by Eq. [T] with 7' = 0. 
Analytic solutions describing a (single) drop contour can 
be parametrized around a local minimum/maximum (s = 



52, x = 0, and ft = h N , see Fig. [H dashed lines) [Iol-[l2|: 
ft^ [Ai(s)Bi(s 2 ) - Ai(s 2 )Bi(s)] 



(3Cb/2)V3p fla ( 5 ) 



ft 7 '=oO) = h N /(7rp S2 (s)) 2 , 



(6) 
(7) 



with p S2 (s) = Ai(s)Bi'(s 2 ) - Ai'(s 2 )Bi(s). For s 2 G 
[—1.01879,0] Eqs. [6] and show a local minimum at 
s = s 2 . The abscissa is scaled by the position of the re- 
ceding contact line x r = h N ( § Ca) _1/ ' 3 Ai(s 2) / Ai' (s 2) , ren- 
dering the plot independent of Ca. The ordinate is scaled 
by ft^ (spherical cap profiles appear distorted). On the 
right, at x = 0, the drop contour is connected through 
a neck of height ft^ to the contour of an infinitely large 
liquid volume: the contour diverges at x > . At the 
branch limit s = s (oo) w -1.732 + 0.595 s 2 (see e.g. [Tl| 
for details), both x 7 /=o — > 00 and ft 7 /=o — > 00. 

We assume now that Eqs. [6] and [7l for x < indi- 
vidually describe the main body contours of drop 2 and 
drop 1, respectively. For composing the twin drop con- 
tour from these two individual contours it is assumed 
further that the contour of drop 2 is modified at x > 
(i.e., to the right of the neck in Fig. |4]) by the local gradi- 
ent 7' (Eq. [5]). In agreement with the experimental and 
numerical results presented above this modification does 
not affect the (gradient-free) contour section, describing 
the main body contour of drop 2 (Eqs. [6] and [3 for x < 0). 
Instead, the local gradient 7' affects the asymptotic cur- 
vature of the contour for x — >• 00, far away from the 
neck region. This gradient-modified curvature is then 
matched to the apex curvature of another (gradient-free) 
drop contour, which describes the shape of drop 1. At the 
apex of drop 1, the curvature is always negative whereas 
the contour of drop 2 from Eqs. I6|7I for x > (without 
gradient) always has a positive curvature. Therefore, to 
match both curvatures, a (local) gradient has to be suf- 
ficiently "strong" to establish a negative curvature for 
x > 0. Then, such matching links the parameters of the 
individual solutions to a twin drop configuration. 

Impact of the surface tension gradient. — The curvature 
at the neck (x = 0) is derived from Eqs. El 



h", =0 (x = 0) = -2(3Ca/2) 2 / 3 52 /ft^ 



(8) 



For x 7 / = o —> 00 (the diverging contour section) the cur- 
vature is finite and positive: 

h% (x -»■ 00) = 27r 2 (3Ch/2) 2 /%( S ^) 2 //i Jv . (9) 

With y 7^ restricted to the neck region at x > 0, the 
curvature change can be approximated: 



00 y 



272 Jo ft 7 '=o 



dx. 



(10) 



Here, the true h(x) is approximated by fty=o, the con- 
tour without gradient. Of course, if we assume a signifi- 
cant impact of the local 7 / (x), at some distance from the 
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FIG. 5: Capillary number for non-coalescing drops with differ- 
ent contact angles (O e =O e ,2=O e ,i) as function of the surface 
tension difference A7/72. 



neck, h will be different from h 7 >=o. However, ^'(x) is 
localized and h(x) grows quickly. Thus, with increasing 
distance from the neck, 7' /hy=o (the integrand) rapidly 
becomes negligible: The integral error remains small, the 
upper integration limit is irrelevant and can be set to in- 
finity [14]. Eq.[10]is integrated analytically by expanding 
hy=o(x) around x = to second order (see Eq. [8]): 



3 r 
W2J0 



-dx i 



7'=0 



(11) 



Here k = Bo(^Ca) 1/s and Ay « A7/72. q(k) is 
lengthy [14], but (for typical parameters) can well be ap- 
proximated in powers of k: q(k) ~ 1 + a/tt^/4 + k(ir + 
fe/3)/16. Eqs. [TO and [TO] yield: 
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. (12) 



The second term indicates that the asymptotic curvature 
may indeed become negative for sufficiently large A7. 

Matching the curvatures. — Eq. [12] estimates the 
asymptotic curvature for x — > 00 as a function of a local 
gradient. As condition for a twin drop configuration this 
curvature shall match the apex curvature of a gradient- 
free drop shape (drop 1). To this end drop 1 is described 
by Eqs. [6] and [7] (renaming S2 to si, si > to parametrize 
at the apex (maximum) instead of the neck). The apex 
curvature is (Eq. [8|) 



til 



-(3Ca/2) 2/3 2s 1 /h 1 , 



(13) 



introducing the apex height hi. 

A moving twin drop configuration with a stationary 
contour means identical capillary numbers Ca in Eqs. [T2l 
and [13] This yields: 

= WCa/2f*q{k) [h^/h, + 7T 2 p' S2 (s^) 2 ] . (14) 

Equation fTTl still contains the parameters 52, si, h N , 
hi (and implicitly Bo). These are linked to experimental 



parameters such as the static equilibrium contact angles 
6 e ,i and Q e ,2, the drop volumes (= areas A\ and A 2 
under the drop contours ) and D. The contour area A2 
(x < 0) is connected to h N by: 

h 2 f°° 

A 2 = (3Ch/2)-V^ / p^syUs, (15) 

Matching the receding contact line to a microscopic 
solution [12|, [l3| gives 52 as a function of Ca, with the 
e> 2 and a slip length A c (besides A2) as parameters: 



Ai'(s 2 ) 
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(16) 



Matching the advancing contact line of drop 1 to a mi- 
croscopic solution (with the inflection point of solutions 
Eqs. [6] and [7] approximated by S{ w s\ — 1) yields [12[: 
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(17) 

Ai for drop 1 following Eq. [15] (with suitable integration 
limits) determines the apex height hi. 

Now the system is closed. The A7 that propels 
the drops with Ca results from inserting solutions to 
Eqs. [[6] and H7] into Eq. [H Figure [5] presents the 
ensuing Ca as function of Ay for various 6 e (assuming 
0e=©e,2=0e,i)- Experimentally this means, given com- 
binations of e , Ay and 77 result in certain Ca i.e., the 
twin drops move at a certain speed. Theoretical (Fig. [5]) 
and experimental (Fig. [2]) findings agree quantitatively 
(within w 10 %). Both show that vn is approximately in- 
dependent from A7 for a wide range of A7 (which is quite 
counterintuitive, A7 is obviously the source for the mo- 
tion). The analysis predicts that the contact angles de- 
termine the order of magnitude for Ca. Assuming exper- 
imentally realistic conditions, other parameters change 
Ca by less than ±25%. As depicted, changing D by a 
factor of 100 (e.g., due to shear-induced dispersion [14]) 
barely affects Ca. Different absolute or relative drop sizes 
also have little influence because the (large) neck curva- 
ture, not the (much smaller) drop curvatures dominates 
capillarity. 

The analysis leading to Fig. [5] assumes stationary con- 
ditions. This requires large drop volumes so that the 
flow between drop 2 and drop 1 changes the composi- 
tion of drop 1 and thus A7 only slowly. Also, A7 has to 
be sufficiently large so that minor changes in Ay barely 
change Ca and vjy is approximately constant (horizontal 
sections of Ca vs Ay in Fig. [5]). The range where Ca 
decreases rapidly with decreasing A7 is quantitatively 
out of the scope of the analytical approach presented 
here. Experimentally, immediate coalescence is observed 
in this region. Nevertheless, via the flow through the 
neck, we can estimate with our approach the lifetime of 
the non-coalescence [14] (which agrees well with the ex- 
periments). 
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Conclusion. — We present new data from the state of 
non-coalescence of two sessile drops from different, com- 
pletely miscible liquids and analyze the hydrodynamics 
of the (temporarily) stationary twin-drop configuration. 
In a thin film approximation the moving twin drop con- 
figuration is described as two moving drops that are 
connected by a liquid neck. Through this neck liquid 
flows from upstream to downstream drop. An advect ion- 
diffusion balance establishes a localized, (temporarily) 
stable surface tension gradient close to the neck, which 
causes a Marangoni flow that sucks liquid out of the neck. 
This counteracts the capillary-driven flow into the neck 
and thus stabilizes non-coalescence. The whole system 
forms a self-stabilizing, traveling wave (twin drop con- 
tour). The theoretical predictions are in quantitative 
agreement with the experimental findings. 
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